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Abstract 

Quantum states evolving, for some < p < oo, into a sequence of mutually orthogonal states 
at times t + nT (n = 0, ...,p - 1) exhibit an interesting physical effect. The analysis of the 



anticipation probabilities p n = 



\Q t +\' qt+nT ) 



2 



shows, that for randomly chosen states q 



measurements of the state at time t + T/2 reveal information about the state at time 
t±nT (n = 0, ...,p - 1), anticipating future states and reflecting past states with significant 

probability. For p < oo and any < S < 1, the probability to measure at time t + T - states from 
the set [q t+n7 -| < \&\ exceeds a 2 S with certainty as p -> oo, where a 2 is a constant 

independent from the period. For p = oo and any N, the probability to measure at time t + T - 

states from the set {q t±nr |n > N} is > a 2 . We characterize the spectrum, establish an analog 
to Planck's relation, define a random sampling scheme, analyze the resulting distribution of 
the anticipation probabilities and point out applications. 
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Introduction 



h . 



Under a Hamiltonian H, a quantum state q evolves into an orbit q t = U t q, where U t = e 
Due to the unitarity of U t , the amplitude d^i q (T) = (q t ,q t +T) = / e~ aT dn q {X) is independent 
from t. If for some integer < p < oo: dfi q (nT) = S nmod p (n e Z), then {q t+nT \n = 0, ...,p - 
1} is a set of mutually orthogonal states, for every t. 

The subject of our study are states with this property, named orthogonal evolution of period p 
and step size T, and the anticipation amplitudes a n = d[i q ((n + 1/2)T), i.e. the inner product 

of states at time nT and T -, as well as the anticipation probabilities p n = \a n \ 2 . It turns out, that 

for randomly chosen states q, measurements of the state at time T - reveal information about 

the states ±nT (n = 0, ...,p - 1). Future states are anticipated and past states reflected with 
significant probability. These probabilities and the strength of this effect, which we name 
Quantum-Mechanical anticipation, are quantified in this paper. 

Chapter 1 shows that orthogonally evolving states have pure point or absolutely continuous 
spectrum dfi q {X), but neither mixed nor singular continuous spectrum. The fundamental 

upper bound hT -1 < -inf Ao (\H\) in terms of the spectral quantity (\H\) is established, which is 

an analog to Planck's relation. The spectral difference is defined in terms of the reduced 
spectrum and used to represent the anticipation amplitudes. 

In chapter 2 we justify from the symmetries of the set of self-adjoint Hamiltonians the random 
sampling scheme underlying our analysis, derive the statistical distribution of the spectral 
difference as q is randomly sampled from all orthogonally evolving states with period p and 
step size T, identify the statistical quantities used in the subsequent chapters to assess 
anticipation strength and show that multiples of p are unlikely to occur by chance. 

In chapter 3 the anticipation effect is analyzed under point spectrum. The distribution of the 
p n and of various observables is derived for some model states as well as under random 

T 

sampling. Theorem 1 states, that for any < 8 < 1, the probability to measure at time - a 
state from the set lq nT n — v -< |<s) exceeds o 2 8 with certainty as p -* oo, where a 2 is a 



constant independent 
expectation « p r 



rom the period. The operator (n r ) (defined in section 2.3) has 



Chapter 4 determines anticipation under absolutely continuous spectrum. Theorem 2 states 
the constant lower bound a 2 for the probability to measure at time - a state from the set 

{q ±nT \n > N}, and infinite expectation for the operator (n r >. The statistics in the continuous 
case are the limiting values of those in the periodic case as p -> oo. 

In the final chapter we draw conclusions, discuss some open questions and point out 
applications. Anticipation in general type quantum evolution will be analyzed in a future 
paper. 
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1 Orthogonal Evolution 

After introducing some terminology, we characterize the spectrum of orthogonally evolving 
states, derive necessary and sufficient conditions for their existence, establish a fundamental 
bound on the frequency and represent the anticipation amplitudes in terms of the reduced 
spectrum. 

Definition 1 

Let H be a Hamiltonian, q a state and < p < oo an integer. Throughout the paper we will 
assume that \\q\\l = 1. 

a) q is evolving under// into q t = U t q, where U t = e a , has spectral measure dn q (X). 

b) q has orthogonal evolution of period p and step size T, iff 

dfiqinT) = S nmod p (n G 1). 

Here dfi q (T) = (q t ,q t +T) = J e~ aT/h dfi q (A) is the Fourier transform of dpi q (X). 

c) a n = dpL q {(n + 1/2)T) is the anticipation amplitude of q nT . 

d) p n = \<x n | 2 is the anticipation probability of q„ r . 

e) The reduction ofdpL q (X) modulo m is dv qm (jc) = Z K= «i mod dfi q (X). 

h 

d[i„{nm 

f) The K-cut of dfi q (A) modulo m is the function c K>m (n) = — - — (n £ 1). 

g) The standard scale is defined to make 7=1. 

1.1 Spectrum 

Orthogonal evolution constrains the spectrum to either pure point or absolutely continuous, 
as stated by the following proposition. 

Proposition 1: Characterization of the Spectrum (Orthogonal Evolution Constraint) 
Let q have orthogonal evolution with period p and step size T. In standard scale, 

a) If p < 00, then v q%ln ({— )) = - (k = 0, ... ,p - 1) is pure point. 

b) If p = 00, then dv q 2n {K) = — [0 < k < 2n] is absolutely continuous. 

Proof: 

-—t 

The spectrum of H and U T are related by U t = e n . U T is a bounded operator with spectrum 
on the complex unit circle [RS80]. Its spectral measure is uniquely determined by 
(f(U T )q ,q ) (/ £ C°°(C)). These functions in turn are uniquely determined by the monomials 
z k , which due to mutual orthogonality satisfy (Uj:q ,q ) = (q nT ,q ) = S nmod p . Thus the 
spectral measure of U T is equally distributed, in the periodic case over the unit roots, in the 
infinite case over the whole complex unit circle. Singular continuous spectrum is excluded. 
Mixed spectrum cannot occur, as in the periodic case the initial state must recur, whereas in 
the aperiodic case dn q (t) e o(l) (see [YL96]). 

The numerical values (in standard scale) follow from $l % e~ inK dv qi2n (K) = 5 nmod p (see 
definition 1b) and the normalization of q. The integral becomes p _1 YJkZl e ~ 2nink /p ar| d 
Q n e~ inK ^| for p < 00 and p = 00, respectively, which is easily verified to yield the correct 
values. 

q.e.d. 
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A more elaborate proof facing the various intricacies and the explicit determination of the 
Fourier expansion of a continuous function / on [Q,2n] such that dv q2n (.K) = D 2 f will be given 
in a future paper. 

Proposition 1 has obvious consequences for the existence of orthogonal evolution. 
Corollary 1 : Existence 

Necessary and sufficient for the existence of a state q evolving orthogonally with period p 
under Hamiltonian H is that either p < oo and the reduced point spectrum of H contains p 
equidistant points (such as the X 2 operator), or p = oo and H has non-empty absolutely 
continuous spectrum (such as the —A operator). 

Proof: 

In either case states meeting proposition 1 are easily constructed. 
1.2 Fundamental Frequency Bound 

In terms of the spectral quantity (\H\) we derive here the fundamental frequency bound 
fij-i < ^(\H\), an analog to Planck's E = hv. 

Definition 2 

For any Hamiltonian H and state q, (\H\) « f\A - A \dn q (A). 
Proposition 2 

Let q be any state evolving under Hamiltonian H. Then 

a) (\H\) < oo iff \(H)\ < oo. 

b) (\H\) > iff q * and q not an eigenstates with zero eigenvalue. 

C) \(q,q t )\ > \Re(q,q t )\ >l-(\H\)t/k. 

d) If dfi^T) = 0, then T > sup Ao (h/(\H-A \)). 

e) The supremum is attained when A equals the median of pi, i.e. At(U ,co]) = 

M([-oo,Ao]). 

f) mfx Q {\H - X \) = iff q is an eigenstate. 

Proof: 

a) By elementary measure theory, (H) = J Adfi q (A) exists iff J 00 Adfi q (A) < co and 

f_ m Adpi q {A)< oo, i.e. iff <|//|> < oo. 

b) As (|//|) is a sum of two non-negative quantities, it vanishes only, if dn q {A) = or 
supp [i q (A) = {0} 

c) By the Schrodinger equation and the triangle inequality, 

||lk t = |Hf q t ,q) + (l.-lfrt)\ = |^ / 1 sin f (A) | < | <|//|>. 

Therefore \\q t - q\\j < \{\H\)t. The result follows from Re(q,q t ) = 1 - \\q t - q\\\/2 

d) By orthogonality, ||q r -q||2 = 2. As T is the same for H - A , the claimed result 
follows. 
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e) From g- (\H - A 1) = - J °° dn q (A + A ) + J °° d H (A - A) = 0. 

f) Elementary. 

q.e.d. 

2c) yields a lower bound on the modulus of the anticipation amplitude at t = T - . shows that at 

given (\H\) the frequency of the orthogonal evolution cannot be arbitrarily high, or 
equivalently the step size arbitrarily small. A bound for the passage time based on (H 2 ) has 
been established by Fleming [FI73] and discussed by various authors. The bound given in 
2d) for the passage time can be strengthened to a bound for the step size of orthogonal 
evolution: 

Corollary 2: Fundamental Frequency Bound 

The frequency of an orthogonal evolution is bounded by 

hT- 1 <U\H\). 

IT 

Proof: Proposition 2d) and 2e) imply hT' 1 < (\H - A |). We conclude from proposition 1 that 
among all Hamiltonians providing orthogonal evolution those with spectrum evenly 
distributed in [-7r,7r] minimize (\H\), which is easily found to be |. 

q.e.d. 



1 .3 Representation of Anticipation Amplitudes 

From now on all spectra are in standard scale. In proposition 1 we have determined the 
reduced spectrum dv q2n {K). Using a notation for the reduced spectrum on the doubled 
interval [0,47r] we establish a representation of a n which we use throughout the paper. 

Definition 3: Spectral difference 

a) Let p < co. Then x kiS = v qAn (^2it [s + j) (5 = 0, 1; k = 0, ... , p - 1). 
7k = x k,o ~ x k,i is tne spectral difference. 

b) Let p = 00. Then x KS = —v qAn (2ns + k) (s = 0,1; < k < 2n). 
y K = x k,o - x k,i is tne spectral difference. 

Proposition 3: Representation of anticipation amplitudes 

a) Ifp < co, then < x kiS < i, ~<y k x k,o + x k,i = J- 

b) Ifp = co, thenO<x K , s <^,-^<y K <^, x Kfi + x Kil = ± 

c) If p < 00, then a n = I^ljy fc e ~ 2ni ^ n ~^ k/p (n integer) . 

d) If p = 00, then a n = j 2r ' y K e~ l ^ n ~^ K ^ (n integer). 

Proof: 
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a) From proposition 1a) and definition 3a). 

b) From proposition 1b) and definition 3b). 

c) By definition 1c) and 3a), 

p— i 

Z- 2 ,fn4W P i —2ni(n—L\l+k/p) 
*fc,o e v 2) i" +x kil e V 2)^ 

k=0 
p-1 

^(xk.o-Xk.iK 271 ^-^. 

k=0 

The claim follows from definition 3a). 

d) By definition 1c) and 3b), 

f 2n -i( n -\) K dK -\\2n+K) dK 

«n=J X Ki0 e V 2) — + XKie V 2)" > — 

[ ln ( \ -i(n-iy dK 

The claim follows from definition 3b). 

q.e.d. 
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2 Random Sampling 

Corollary 1 implies that orthogonally evolving states exist under infinitely many dynamics. In 
the remaining chapters we investigate, how likely such states are to exhibit anticipation of 
significant strength. To this end, we define in this chapter a random sampling scheme which 
will underlie all our subsequent analyses. We justify this scheme by the symmetry properties 
of the set of self-adjoint Hamiltonians. Under this scheme we derive the statistical distribution 
of the spectral difference, which determines the anticipation amplitudes. We define the set of 
statistical quantities used in the subsequent chapters to assess anticipation strength, and 
finally show that multiples of the period are unlikely to occur by chance. 

2.1 Sampling Scheme 

In order to investigate the likelihood of physical states to exhibit anticipation of a certain 
strength we must consider the set of all orthogonally evolving states and the way such states 
arise. We focus here on a laboratory situation, where the experimentalist prepares an 
experimental setup and an orthogonally evolving target, and then studies its evolution under 
the laws of nature. The setup is modeled by a Hamiltonian and the target by a physical state. 

What Hamiltonians and states are preparable depends on the experimentalist's skills and 
capabilities and the constituents available in nature. Abstracting from these practical 
limitations, we only constrain the experimentalist to prepare self-adjoint Hamiltonians and 
states exhibiting orthogonal evolution, and analyze the outcome if these are randomly 
chosen, letting freely fluctuate all degrees of freedom except those fixed by the orthogonal 
evolution constraint (proposition 1). 

We are going to define a sampling scheme, where, for a given step size and period, first a 
self-adjoint Hamiltonian is chosen from the set of all self-adjoint Hamiltonians admitting 
orthogonal evolution, and then a state is chosen from the set of states orthogonally evolving 
under that Hamiltonian. 

Proposition 3 shows that orthogonal evolution and anticipation are completely determined by 
the spectral properties of the Hamiltonians and states involved, more specifically by the 
spectral difference. Thus, for our purpose, only constraints on and degrees of freedom of 
spectra need consideration. 

Inverse theory tells us that, for any bounded real measure fi and any positive integer d, there 
are uncountably many self-adjoint operators H and states q e L 2 (U. d ) with spectral measure 
ju. This holds as well for at least one Schrodinger operator -A + V with self-adjoint extension 
H. This enables us to define the sampling scheme in terms of probability spaces of spectral 
measures. 

The core feature of spectral measures of orthogonally evolving states are the K-cuts modulo 
2n (see definition 1f), related to the spectral measure \i by dfi(2nn + k) = c K2n (n)dv q2n (K). 
The quantities x ks from definition 3 are x ks = Y 1 c K2n {2n + s)dv q2n (jc). 

Whereas the reduced spectral measure is fixed by the orthogonal evolution constraint, the k- 
cuts are arbitrary doubly-infinite sequences with non-negative values, summing up to 1. 

Let's firstly consider their supports suppc K2n , i.e. the patterns of points neZ where 
c K,2n( n ) ^ 0- This pattern is a property of the Hamiltonian, whereas the size of c Ki2n (n) is a 
determined by the state. 
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Self-adjointness allows for any such pattern at any index k, and to choose patterns at 
different indices k independently. This holds for point spectrum as well as for continuous 
spectrum, with the only constraint that for fixed n the dependence from k must be absolutely 
continuous. While the following arguments will be given for point spectrum, they are easily 
extended to continuous spectrum as the limiting case. 

This symmetry of the set of self-adjoint Hamiltonians allows us to require the sampling 
scheme to make the supp c K2n a family of independent identically distributed random 
variables. As these symmetries hold as well when constraining the Hamiltonians to admit 
orthogonal evolution, we demand the same property in this case. 

Given any Hamiltonian, the value of c Ki2n (n) on suppc K2 „ depends only on the state q. 
When randomly sampling states q e L 2 (R d ) meeting the orthogonal evolution constraint, 
then c K 2n can attain any non-negative function f e I 1 (supp c K 2n ) meeting this constraint, 
independently for any index k, because orthogonal components of these states belonging to 
different indices can be sampled separately in the discrete case. As by the preceding 
requirement the supports are i.i.d., we impose the same requirement on the sampling 
distribution of the K-cuts. 

As all distributions as well as the reduced spectrum are index-independent, we will omit 
indices in the following. Let 6 p be sampling distribution of the K-cuts c K 2n at period p. 

Notice that 6 p is not obtained by conditioning from another, unconditional distribution. The 
orthogonal evolution constraint is met by normalization of the K-cuts, not by conditioning of 
some distribution. Practical sampling procedures will be designed to meet the constraint 
approximately, but not by selection of "good" states from unconstraint samples. Only the 
latter procedure would yield a conditional distribution, whereas the real procedure is 
sufficiently modeled by a procedure where orthogonal components are sampled and 
normalized independently (index by index). 

The period p only determines the orthogonal evolution constraint. As this constraint is met by 
normalization, it only enters into 6 p as a scaling factor: 9 p (c Ki2n ) = q (^c K2n ^ for any p, q 

and c K 2n . As our reasoning did not depend otherwise from the period p, we finally postulate 
the sampling distribution of fc-cuts to be independent from p, up to the orthogonal evolution 
constraint and the scaling factor. 

Our requirements are summarized in the following definition. 
Definition 4: Random Samples 

A family of sets S p of measures a is a family of random samples of measures of orthogonally 
evolving states with period p < °o iff the following holds for all p: 

a) The reduced spectrum v q 2n meets proposition 1 above. 

b) For p < co, the K-cut-valued random variables c K are i.i.d. under probability measure 
6 p (c) = 6 l (pc) independent from k. 

c) For p = oo, the K-cut-valued random variables c K are i.i.d. under probability measure 
0oo(c) = 0i(27rc) independent from k. 
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From the preceding arguments and definition 4 follows: 
Corollary 3: Existence 

The class of self-adjoint Hamiltonians admits for sampling schemes yielding random samples 
as defined above. 

Further evidence justifying our sampling scheme will be provided by the main results and 
discussed in the final chapter. 

2.2 Distribution of the Spectral Difference 

Random sampling according to the above scheme yields a probability distribution of the 
spectral difference function y k . It is no surprise that these are again i.i.d.. 

Proposition 4 

Consider the sampling distribution as q is randomly sampled under the sampling scheme of 
definition 4. Then for some distribution function F(z) on f-1,+1] the following holds: 

a) If p < oo, then y k (fe = 0, ...,p — 1) is a sequence of i.i.d. random variables with 
distribution F(py) with density p/(py) for (-^ < y < ^. 

b) If p = oo, then y K (0 < k < 2u) is a family of i.i.d. random variables with distribution 
F(2ny) with density 2nf{2ny) for (-^ < y < 

c) If the random sampling scheme is invariant under translations of the energy scale, 
then the density / is an even function. 

Proof: 

a) For p < oo the distribution of x ks is X s (x) = 0i({c 2jr |p:x: > Zc 2;r (2n + s )})> f rom which 
F is obtained. The scaling is due to proposition 1. The i.i.d. and period dependence 
are direct consequences of definition 4. 

b) For p = oo the distribution of x kiS is X s (x) = Oxtfan^nx > X c ln {2n + s)}), from which 
F is obtained. The scaling is due to proposition 1. The i.i.d. and period dependence 
are direct consequences of definition 4. 

c) Translations of the energy scale by 2nm effect a shift by —m on the K-cuts, 
transforming them to c K 2n (n - m). If m is odd, the sign of the spectral difference is 
toggled. 

q.e.d. 



2.3 Anticipation Statistics 

We identify here the statistical quantities used to assess anticipation strength in the 
subsequent chapters, and recall some statistical terminology and basic facts. 

The following anticipation statistics will be analyzed in the subsequent chapters, with n as by 
definitions 5a and 5b below. 

a) p n , the probability to measure q nT at time T -, and p N = P(n > N). 
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b) Ptot = Z Pn . tn e probability to measure one of the states q nT at time T -. 

c) Observable n r (r > 0), with value n r in state q nT (n > 0) and expected value 

We perform the analysis for minimum and maximum anticipation model states as well as 
under biased and unbiased random sampling according to definition 4 above. In the sampling 
scenarios, we provide the sample mean and variance as well as stochastic asymptotics of 
the probabilities. 

The following definitions will be used: 



Definition 5 



a) If p < oo then n = f N m ° d P (N m ° d P ~ (P + 1)/2) 

{p + 1 — (|n| mod p) (otherwise) 

b) If p = oo, thenn = \n\. 

c) m n = E(z n ) (n = 0, 1, ... ) is the n th moment. 

d) o 2 (z) = E((z — mi) 2 ) is the variance. 

e) z E0 p (1) iff Ve > 3c: lim sup P(\z\ > c) < s. 

f) z £ o p (l) iff V£ > 0:limP(|z| > e) = . 



p and o p are Pratt's stochastic asymptotics [Pr59], the former expressing convergence in 
distribution, or stochastic boundedness, the latter stochastic convergence to zero [BFH75]. 

Recall the following facts: 



Proposition 5 

a) If the probability density is an even function, then m 2n+1 = 0. 

b) E((cz) n ) = c n m n . 

c) a 2 = m 2 — m 2 . 

d) P(|z|>(£(|zn) 1/? A)<^ (r>0,A>0). 

e) P (^-^ > a) < where z = m 1 . 

f) z = z + O p {o). 

g) If a £ o(z) then | = 1 + o p (1). 



Proof: 

Parts a-c) are elementary, d) is the Bienayme-Chebychev inequality [TJ71], e) its special 
case r = 2, f) and g) direct consequences of e) and definition 5e-f). 



2.4 Multiple Periods 

Random sampling of orthogonally evolving states with period p does admit states with period 
np for some n > 1. Such states have y k = V/c, thus = 0. However, they are not likely to 
occur by chance in random samples. 
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Proposition 6 

Let q be a state with orthogonal evolution of period p. 



a) p < co. For any s > 0, let n be the cardinality of the set {k\ \py k | < s, < k < p}. Then 
n has the Binomial distribution with parameter p and F(e). 

b) p = oo. For any s > 0, then the set {?<:||y K | < s, < k < 2n) has Lebesgue measure 
F([- £ ,+ £ ]). 

Proof: 

a) Immediately from proposition 4. 

b) From proposition 4, as the variables are i.i.d. 

q.e.d. 

This result rules out the occurrence by chance of period np or, equivalently, step size T/n 
when randomly sampling states with period p and step size T. However it does not guarantee 
strong anticipation. This will be provided by the proof of constant probabilistic lower bounds 
on p N in the next two chapters. 
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3 Anticipation (Point Spectrum) 

Here we determine, for model and random states with p < oo, the values and distributions of 
the anticipation statistics identified and defined in section 2.3, and prove theorem 1. 



3.1 Minimum Anticipation 

We consider here the extremal case of constant y k . Notice that these model states have zero 
occurrence probability. 



Proposition 7 

Let py k = y (k = 0, p - 1) for some -1 < y < 1. Then 

a) p n = —, — n— with maximum value ~ at n = and n = 1, minimum value 

p 2 sin 2 n(n— jJ/P 71 

« p~ 2 y 2 at n = ||j and p„ £ 0(n~ 2 ) (n -> 

b) p N 6 O^- 1 ) (ff f||). 

c ) Ptot = y 2 - 

d) (n 1 ) = lnp + 0(1), (n r > = 0(p r - 1 ) (r > 1). 



Proof: 



a) By proposition 3c) and the premise 

p-i 

-2m (n— jjk/p 



k=0 

V y = 1 — W, = { y 7 T\ — (2n - 1 * p) 

e -2ni{n-f)/p _ 1 e -2m{n-2)/p _ ± p s i n „ (n - ±J/p 

p _1 y (otherwise) 
The maxima, minima and asymptotics are elementary. 

b) From a) by the Euler summation formula. 

c) From Lemma 1 below. 

d) From a) by the Euler summation formula. 

q.e.d. 



Lemma 1 

p— i 

p~ 2 sin -2 n(n — ^j/p = 1 

Proof: 

Let q be a superposition of exactly p eigenfunctions. Then Span({q t \t e E}) = Span({q kT \k £ 
Z}). Therefore p tot = \\q T /2 \\ 2 = 1- As in this case \y \ = \py k \ = 1, the claimed result follows. 
The alternatively, direct summation, is easily achieved by adapting theorem 4.9a of [He88] to 
periodic functions. 
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3.2 Maximum Anticipation 

We consider here the extremal case of alternating y k = {-l) k y. Notice that these model 
states have zero occurrence probability. 

Proposition 8 

Let p be even, py k = {-l) k y (k = 0, ... , p - 1) for some -1 < y < 1. Then 

a) p n = — jt — , which has maximum value ~%- at n = [^1, minimum value 

p 2 cos 2 n[n--)/p n 12 1 

« p _2 y 2 at n = and n = 1, and p n £ 0(p~ 2 ) (n £ o(p)). 

b) Pn = Ptot ~ 0(P _1 ) (p e o(p)). 

c) p tot = y 2 - 

d) (n r > = 0(p r +p r_1 In p). 

For odd p, this evolution degenerates to an orthogonal evolution with period p and step size 
T/2. 

Proof: 

a) By proposition 3c) and the premise 

p— i 

a n =p-i y ^(-l)fc e - 2 -("4>/P 

k=0 

( -1 l-t-i)>e- 2 < n -l) = 2p-V t 

= \ i+e _ 2ff i(„-l)/p 1+e -2,ri(n-i)/p 

lo 

The maxima, minima and asymptotics are elementary. 

For odd p and y = ±1, the support of dv^O) consists of p equidistant points. Thus 
the evolution is of minimum anticipation type with period p and step size T/2. 

b) From a) by the Euler summation formula. 

c) By Lemma 1. 

d) From a) and the Binomial theorem one obtains (n r ) ~ £y, n (|) ' ( _n ) 7_2 y7. f rom 

which the result follows by the Euler summation formula. The notation a- denotes the 
b th falling power of a. 

q.e.d. 

3.3 General Case 

We consider now sampling distributions such that y k {k = 0, ...,p - 1) is a sequence of i.i.d. 
random variables with distribution F p (y) (—^<y<^j with moments p~ n m n . 



e -ni(n-z)/p 

= y , p- (p even) 

p cos n\n—-\/p 

(p odd) 
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Proposition 9 



a) E(p n ) = p- 1 a 2 + 



2 

mf 



p 2 sin 2 n(n—^j/p 



= p~ l a 2 + f 



0(p" 2 ) (n £ n(p)) 
0(fT 2 ) (n G o(p))' 



ECPn) has maximum value p 1 c 2 + — m 2 at n = and n = 1. 
^CPn) has minimum value p _1 o- 2 + at n = |"||. 



b) p n = £(p n ) + - 



(= K = m? ) 

eO p (p~ 1/2 ) (2n-l = p) 

G O p (p" 1/2 n _1/2 ) (n£o(p)) 

£ p (p _1 ) (n £ H(p), 2n - 1 * p) 

leOpCp- 1 ) ( mi 6 0(p- 1/2 )) 



c) E(p tot ) = m 2 . 



d) E(p ]V ) = (l-2p- 1 iV)a 2 + m 2 ^ 



(JV = 0) 

(tf G o(p),iV * 0). 



^p — 2p 1 JV) (otherwise) 



= 



e) p w - £ (P^ + ( e0p ( p -i/2) (ot/iervv ; 5e y 



f) £■ (n r ) = + OCp^- 1 ). 



Proof: 

a) From proposition 3c) and the premises, as detailed in the appendix. 

b) From proposition 5f, where Var(p n ) is calculated in the appendix. 

c) See appendix. 

d) See appendix. 

e) From proposition 5f, where Var(p N ) is calculated in the appendix. 

f) From proposition 8a and the Euler summation formula. 

q.e.d. 



Proposition 9 shows that anticipation is fairly strong. Whereas the expected probability is 
independent of the period, the standard deviation is 0(p _1/2 ). We summarize these findings 
by 
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Theorem 1: Anticipation (Point Spectrum) 

a) When randomly sampling with mean m 1 and a 2 , then, for any < 5 < 1, 
measurements at time T - yield a state from the set {q nT \N<n<p + l- N} with 
probability p N = a 2 8 + o p (l) for N = |(l - <5). 

b) For any < 5 < cr 2 and any JV, the probability p N > 8 with certainty as p -» oo. 

c) The operator (n r > has expectation 0(p r ). 



Proof: 



a) By proposition 9b. 

b) By a) and definition 5f). 

c) By proposition 9f). 

q.e.d. 

Notice that for unbiased sampling m 1 = and a 2 =m 2 , which by proposition 4c) is easily 
achievable. 
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4 Anticipation (Continuous Spectrum) 

Here we determine, for model and random states with p = oo, the values and distributions of 
the anticipation statistics identified and defined in section 2.3, and prove theorem 2. 

4.1 Minimum Anticipation 

We consider here the extremal case of constant y k . Notice that these model states have zero 
occurrence probability. 

Proposition 10 

Let y K =y (0 < k < 2n) for some -1 < y < 1. Then 

a) Pn = y 2 7r -2 (n — with maximum value ~ at n = and n = 1. 

b) pjv e 0(N~ 1 ) (N ^ oo). 

c) Ptot =y 2 - 

d) (n r > = oo, ( r > 1). 



Proof: 

a) By proposition 3d) and the premise 

f 271 -i( n -^\ K dK -1 / 

a n = I ye v 2; — (n integer) = —iyn 1 1 n — - J 

The maxima, minima and asymptotics are elementary. 

b) From a) by the Euler summation formula. 

c) By a well known result, first discovered by Euler. See also [He88]. 

d) From a). 

q.e.d. 



4.2 Maximum Anticipation 

We consider here the extremal case of y k alternating between y and -y in N equidistant 
intervals. Notice that these model states have zero occurrence probability. 

Proposition 11 

Let y K = (-iy'y < k < 27r ^ +1) , < j < m) for some even < M < 00. Then 



a) p n = 



y tan v 2J 



M 



with minim value ~ y 2 M 2 at n = and n = 1, and maximum 



value ~ 47T 4 y z at n = —. 



b) Ptot = y 2 

r= y 2 (N = 0) 

C) p N Uy 2 (l-0(NM- 2 )) (/VGo(M)) 

(o(^ _1 ) (N -> 00) 

d) (n r > = 00, ( r > 1). 
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Proof: 



a) By proposition 3d) and the premise 



^=Z ( - 1)7y k e 1 j to 





b) Either by direct summation using theorem 4.9a of [He88], or using the fact that the 
functions e~ inK (n £ T) are an orthogonal basis of £ 2 ([0,2n]). 

c) From a) and b). 

d) By comparison with the minimum anticipation case. 



For odd p, this evolution degenerates to an orthogonal evolution with period p and step size 
T/2. 

4.3 General Case 

We consider now sampling distributions such y K (0 < k < 2n) is a family of i.i.d. random 
variables with distribution F(2ny) with moments m n . 

Proposition 12 



a) E(p n ) = 




b) £(Pt t) = ™-2- 



q.e.d. 



c) 




(iV = 0) 
QV ^ oo) ' 



d) p N = E(p N ) with probability 1. 



e) £(n r ) = oo, (r > 1). 



f) These values are the limits of the corresponding values as p -> oo. 
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Proof: 



) E (Pn) = ^J? Q n E{y{K)y{K))e< n -^ K -^dKdK 



2 



u2(n-i) 



b) Ptot - -co 4tt 2 jo Jo 

= lim^oo ^ Q n Q n y{K)y{K)e~ l ~ smm ^ K ~ K ^ dK dK 



where taking the principal value is justified by the fact that p tot < 1. For further 
evaluation, we first notice that the functions h m (x) = sm T (-rt < x < n) have a 

sin- 

removable singularity at the origin and asymptotically vanishing integral for x e 
n(m _1 ), whereas for x £ o(m _1 ) the approximations 

msinmx x 

h m (x) = ynx + = wicosm- + o(l) 

2sin^- ^ 



hold. Approximating cosx in the interval [y-,f 

6.32 from [Ru73] yields a representation 
functional, such that 



by test functions and applying theorem 
\m m _, m mcosm = 2S X of Dirac's delta 



lim ^- ( +n f(x)m cos m^-dx = 2/(0). 

m->oo 2n —ft 2 

The above arguments and the properties of h m imply that 

.k—k 

lim e~ l ~ h m (x) = 5 K _ K > +8 K _ K ' +nmod 2n . 

m-»co 

Therefore 



2n 2n 



Ptot = 4^2 / / + <W+7r mod 2n) dKdK ' 



= — J ^(y 2 ^) + y(K)y(fc + n mod 2tc)) dK, 

from which the result follows by the independence property. 

c) By a) and b) and the Euler summation formula. 

d) Using the same methods as in the proof of b) above with the necessary adaptations 
we obtain 

E (p 2 N ) = 

I 2n 

= E \ Z TeW 3'W3'(Oy(K")y(ic">" i ^) (K " K ' ) " i ( ^ (K " h K d K ' d K " d K '" 



n,o>N 
2n 



= E 



^2 j yWy( K ') ( S k-k ~h M (K- K))dKdK 



= E 2 (p N ). 
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Therefore Var(p N ) = 0. The result follows from proposition 5f. 

e) By comparison with the minimum anticipation case. 

f) By comparison with proposition 9. 

q.e.d. 

There is a seeming contradiction between 9a) and 9b), as Y,Pn = m i ^ m i = Ptot> which is 
due to the fact that lim m £ a„ m does not necessarily equal £ lim a n m . The consequences are 

m 

discussed in the final chapter. 

Proposition 12 proves that anticipation strength in the continuous case is as strong as in the 
discrete case. The p N are of the same size, whereas the standard deviation is zero. 

Theorem 2: Anticipation (Continuous Spectrum) 

a) When randomly sampling with mean m l and standard deviation a, then 
measurements at time - yield a state from the set {q nT \N < \n\ < oo} with probability 

p N =m 2 - ml (l - O^N' 1 )) for all < N < oo. 

b) p N >m 2 — ml with certainty for all N. 

c) The operator (n r ) has infinite expectation. 

Proof: 

a) By proposition 12d. 

b) By a) and definition 5f). 

c) By proposition 12f). 

q.e.d. 



Notice that for unbiased sampling m x = and a 2 =m 2 , which by proposition 4c) is easily 
achievable. 
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5 Conclusions 

5.1 Summary 

We have analyzed the statistical distribution of the anticipation probabilities p n and the 
cumulative probabilities p N when randomly sampling orthogonally evolving states with period 
p. The main findings are that, in the periodic case for any < 5 < a 2 and N = pS, the 
probability p N > S with certainty as p -> oo, and in the aperiodic case p N > a 2 with certainty 
for all N, where a 2 is a constant independent from p. Somewhat surprisingly, the expected 
values of the single probabilities p n have the lower bound p~ x a 2 independent of n, which is 
the main reason for the unexpected strength of anticipation. 

5.2 Statistical Premises 

Our analysis relies on a random sampling scheme justified by symmetries of the set of self- 
adjoint Hamiltonians described in section 2.1. Thus our results hold for random sampling 
from this set of Hamiltonians of which the set of physical Hamiltonians is a sub-set. Lacking 
an exhaustive characterization of the latter, we have based our analysis on the former. 

Our results depend on three features of the distribution of the spectral difference (proposition 
4): Independence from the index k, statistical independence and independence from the 
period p. Let us discuss them one by one. 

Independence from the index follows from the symmetries of the set of self-adjoint 
Hamiltonians. It is hard to imagine what physical effects should break these symmetries. 
Noticing that shifts of the energy scale shift the index k, how should any asymmetry arise. 
We assert that independence from the index holds as well for the set of physical 
Hamiltonians. 

Statistical independence of the y k is provided by separately sampling the orthogonal 
components of a state for different indices k. Though this may not be practical, it shows that 
the assumption of statistical independence is legitimate as a model of reality and 
approximation to experimental practice. 

Independence from the index and statistical independence imply that the anticipation 
strength mainly depends on the variance a 2 , as established in chapters 3 and 4. 

Independence from the period may indeed be broken by some law of nature, or by the 
particular techniques used by the experimentalist. But what will this effect? Will it decrease or 
increase the variance? How and why should a physical or a preparation effect decrease the 
variance of a single spectral difference y k below p~ 2 a 2 as p -> oo, which is the value we have 
assumed? The central limit theorem does not apply in this case, and conditioning has also 
been ruled out. 

As precision is rather expected to decrease than to increase, we maintain our model as a 
good approximation to physical reality and assert that our findings are valid for physical 
Hamiltonians and practical realizations. 

5.3 Features and Limitations 

Orthogonally evolving states are strange objects, as they seem to flip back and forth between 
present and distant times every T - time. The observable <(n - t) r ) oscillates between at full 

steps and 0(p r ) at half steps. Measurements at T - may yield states at distant times. 
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However, this behavior is just due to coherence, not to any non-locality in space or time. 
Future measurements do not affect the state at time T -, nor does measurement at that time 
affect the past state. 

The oscillations may effect translations in space and energy transport. Their size is thus 
subject to energetic and relativistic limitations. 

As an example, periodic motion under a Schrodinger Hamiltonian is accompanied by an 
oscillation of the kinetic and potential energy. Thus measurements of (-A) at full and half 
steps will exhibit such oscillations, their amplitude being related to T - [V, —A] by the mean 
value theorem of calculus. 

In the continuous case, (n r > is related to {X r ) n , as states with absolutely continuous 
spectrum are known to spread to infinity [KL99, KKL01,Si90]. Therefore (X r ) is infinite for 
orthogonally evolving states. <(n - t) r > is not measureable at time T -, because the state is not 

localized, whereas measurement equipment is. However, successive measurements of their 
presence at full and half steps will reveal a difference of relative size a 1 , thus proving the 
disappearance of the target. 

5.4 Applications 

Preparation of orthogonally evolving states followed by measurement after one half-step is a 
generic method for rapid state preparation. 

This obviously applies to quantum computing, which is nothing but preparation of the result 
state from the initial state. It has been explained in [TH0801] that quantum computers 
simulating computations of deterministic Turing machines evolve orthogonally, and how the 
anticipation effect can be exploited to obtain the results of arbitrarily long computations at 

time \ 

5.5 Open Problems 

To our knowledge, orthogonal evolution and the associated anticipation effect have so far not 
been a topic of physical literature and experimental practice. One reason may be, that this 
effect is neither of interest nor easily observable in fast processes of short periods, which are 
predominant in microscopic systems. 

This is different for quantum computers. Their first generation is likely to be of hybrid type, 
with a quantum memory and a classical finite state control and thus not exhibit anticipation. 
Anticipation will be achieved once not only the state but also the finite state control are in a 
coherent state, thus implementing a time-invariant Hamiltonian. We will outline this topic and 
its applications to hypercomputing in future papers. 

The first step towards verification of anticipation will be to prepare quantum systems with 
long periods, and to prepare them such that they evolve orthogonally. Composite systems of 
p quantum oscillators with suitable frequencies are candidate realizations. Once this has 
been done, the next task is the measurement of the predicted anticipation effect. 

Possibly, nature is ahead of us once again, and is using anticipation to speed-up processes. 
It would thus be worthwhile to seek for physical phenomena arising due to anticipation. 
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An obvious question is whether anticipation is specific for orthogonal evolution. We will 
analyze anticipation in general type quantum evolution in a future paper demonstrating that 
anticipation is a universal phenomenon. 
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7 Appendix: Proofs of Proposition 8 
7.1 Mean and Variance of p n 

The tilde operator is defined in definition 5a-b. Define 

2N 

u n = 1 - — 
n' n = 1 - np' 1 

s k M s k (n) <= f p-ie- 2 < n -$ k/p ,s = \s k \= p-K 

S « S n *? ) s fc = p- 1 

_ ie ^( n 4)/p f« 4r (n £ o(p)) 



P sinjr(n4)/p | £ (p-i) (otherwise) 

which satisfies 



■^n "-*p+l— n 



and by Lemma 1 

Z,|5| 2 = i. 

Using 

'"-P e -2n m /p_ 1 (o ( 



*! „-2 e 2 " m ~ 1 _ (P 1 (n = mod p) 

t/ierwise) 

r = |T 2? i_i | which for even p is always zero, 



'n+o-l 



we verify easily 

2s ; (n)s fc (o) = S„S , 

Es;(n)s_ k (o) = S„S 
^ S/ (n)s ; - (o) = T, 
^ Sj (n)s_j (o) = 7;_ 
^ s y (n)s fc (o) = S„S - r n+0 _!, 

Ey** 5; (n)s_ fe (o) = S n S - 7;_ . 
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Now by proposition 3c) 

p— 1 



p-1 - \ 

£ y fce -^("4>/P = p -2 ^ (m2 _ m 2) + £ m 2 s ._ A 
fc=o / fc y,fc 



fc=o 

= p~ 1 (m 2 — ml) + mlS 2 . 
a 2 =m 2 - ml is the variance and thus non-negative. 

E(p tot ) = E(Z n Pn) = T. n E(Pn) and Lemma 1 immediately yield 
E(Ptot) = ™ 2 . 
Again by proposition 3c) 



E(p 2 ) = E 



which we evaluate by partitioning the index set and collecting the contributions of each sub- 
set: 



a) j — k — I — m p~ 3 m A 

b) j k = I = m p~ 2 {S 2 — p _1 )m 1 m 3 
k j = I = m 

I ^ k = j = m 
m k = I = j 

c) j = k I = m p~ A p-m 2 
j = m k = I 

d) j = I k = m (T 2 — p~ 3 )m 2 

e) j = k I ^ m (p _1 — 2p _2 )(S 2 — p~ 1 )m 2 m 2 , 

using {/ = k ^ l,m} — {j = k ^ I = m} 



where the first bracket yields p 2 £(s - s ; )(5 - s,) 
together p _2 (p - 2)S 2 + 2p -3 - p -2 

j ^ k ^ l = m using {I = m ^ j,k} - {I = m ^ j = k} 

j = m=t k =t I using {/ = m =t k, 1} - {j = m =t k = 1} 

k = I ^ j ^ m using {k = I ^ j,m} — {k = I = m} 

f) j = l±k±m ((T — 2p _2 )S 2 + 2p -3 - T 2 )mlm 2 , 

using {j = I ^ k,m] - {j = I ^ k = m} 
where the first bracket yields 

£ sf (5 - 5_ ; ) 2 = ]T s, 2 (S 2 - 25s_ ; + sij ) 

together (T - 2p~ 2 )S 2 + 2p~ 3 - T 2 
k = m j I using {k = m ^ j,l} — {k = m ^ j = 1} 

g) j±k±l±m (S 4 + (8p- 2 - 27 - 4p _1 )S 2 + T 2 — 6p~ 3 + 2p- 2 )m\, 

using {j=tk,l,m}-{j=tk=tl = m}-{j=tk = l=tm} 
{j ± k = m* l}-{j * k = l = m} 

where the first bracket yields £ sj (S - sj ) (S - ) 

= ^ s j s3 — { s f + 2p" 2 )S 2 + 3p- 2 5 y S - p- 4 
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= S 4 — (T + 2p~ l - 3p~ 2 )S 2 - p~ 3 
Together 

S 4 — (7 + 2p~ l - 3p~ 2 )S 2 - p~ 3 - 2(p -1 - 2p~ 2 )(S 2 - p _1 ) - 

((7 - 2p~ 2 )S 2 + 2p -3 - T 2 ) - p~ 2 (S 2 - p _1 ) 

= S 4 - (27 + Ap- 1 - 8p- 2 )S 2 + T 2 — 6p- 3 + 2p~ 2 



This yields 



£(Pn) = V 3m 4 + 4p 2 {S 2 — p 1 )m 1 m 3 + 2p 4 p-m 2 + (7 2 — p 3 )m| 
+ 4(p _1 - 2p _2 )(S 2 - p~ 1 )m 2 m 2 
+2((T - 2p~ 2 )S 2 + 2p- 3 - 7 2 )m 2 m 2 

+ (S 4 + (8p -2 - 27 - 4p~ 1 )S 2 + T 2 — 6p~ 3 + 2p~ 2 )m 4 
> -3 m 4 + 4p _2 (S 2 — p _1 )m 1 m 3 + (T 2 + 2p~ 2 — 3p~ 3 )m 2 
-((27 + 4p _1 - 12p _2 )S 2 + 12p -3 - 4p -2 - 2T 2 )mlm 2 

+ (S 4 + (8p- 2 - 27 - 4p _1 )S 2 + 7 2 - 6p- 3 + 2p~ 2 )m 4 



_ 3 

= p - ■ 



The variance is 



Var(p n ) = E(p 2 ) — p~ 2 (m 2 — m 2 ) 2 — m 4 S 4 — 2p~ 1 m 2 (m 2 — m 2 )S^ 

= E(p 2 ) — p~ 2 m 2 + 2(p -2 — p~ 1 S 2 )m 2 m 2 + (2p _1 S 2 — p -2 — S 4 )m 4 
= p~ 3 m 4 + 4p _2 (S 2 - p _1 )m 1 m 3 + (7 2 + p~ 2 — 3p~ 3 )m 2 

r _ i . _ _i . _ i _ n i\ n 



in 4 t t F w - p 1 )m 1 m 3 + (7 Z + p" 
+ ((27 + 2p _1 - 12p _2 )S 2 + 12p -3 - 2p" 



3p 3 )m| 
2 - 27 2 )m 2 m 2 



+ ((8p- 2 - 27 - 2p _1 )5 2 + 7 2 - 6p- 3 + p~ 2 )m\ 



One easily verifies 



Kar(p n ) 



r= 

£ 0(p _1 ) 
G 0(p _1 n _1 ) 
£ 0(p- 2 ) 

le 0(p- 2 ) 



(m n = m") 
(2n — 1 = p ) 
(n £ o(p)) 

(n £ H(p), 2n - 1 * p) 
(mx £ 0(p- 1/2 )) 
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7.2 Mean and Variance of p N 

We define 

n>N [p 1 n N (otherwise) 

and notice 

From the preceding results we find 

E(Pn) = ^ E (Pn~) = n N (m 2 - ml) + mlU N 

\ (N = 0) 



n>N 



n N o 2 + ml 



(N £ o(p),N * 0) 



\Tl 2 N 

[p~ 1 'n: N (otherwise) 



\\n>jV / / n,d>N 

= £"(Zji, ->iv Zy,fc,i,m yy yfcyiymSy 00 s_ k (n)s l (o)s_ m (o)), 



which we evaluate as before by partitioning the index set and collecting the contributions of 
each sub-set: 



a) j = k = I = m p -1 it^m 4 

b) j k = I = m p _1 7T JV ([/ JV — n N )m 1 m 3 
k j = I = m 

I k = j = m 
m k = I = j 

c) j = k ± I = m n-^iTfjml 

d) j = m ± k = I p _1 7T JV (l — n N )m\ 

which follows from 

Y\ T n-o\ 2 -P" 3 = P~ l n N -p~ l n 2 N 

e) j = I k = m p _1 7r JV (l — n N )m\ 

f) j = k =t I ± m n 2 n N (U N — n N )mlm 2 

using {j = k + l,m] - {j = k =t I = m} 
where the first bracket yields 

P~ 2 Z (S - sj (o)) (S~ - s_j (o))=(p- 1 - 2p- 2 ) |S 2 1 + p- 3 
summing up to 
n' 2 n N U N + p~ x n\ 

j k =t l = m using {I = m =t j,k} - {I = m =t j = k} 

g) j = m^k=tl P _1 (l - 27r w )((/ w - n N )mlm 2 
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using {j = m^k,l}-{j = m^k = l} 
where the first bracket yields 

^ Sj (,n)s_j (o) (s n - s_j (n)) (s o - sy (o)) 

= T n _ S^S - p~ 2 \S n | 2 - p- 2 \S 1 2 + p- 3 
which sums up to 

p _1 f/jv - 2p~ 1 n N U N + p~ l nl = p" 1 y JV (l - 2n N ) + p~ x n\ 
k = I =t j =t m using {k = I =t j ,m} - {k = I =t j = m] 

h) j = I k ± m p _1 (l — 27T JV )(f/ JV — 7T JV )m 2 m 2/ 

using fj = I =t k,m} - {j = I =t k = m} 
where the first bracket yields 

^ Sj (n)sj (o) (s n - s_j (n)) (s o - s_j (o)) 

= £ rn+o-l^o-p-'l^l-p-'l^l+p- 3 

n,o>]V 

summing up to 

p _1 (f/iv - 2n N U N ) + p~ l nfi 

k = m j I using {k = m ^ — {k = m ^ j = 1} 

i) j^k^l^m (ujj + U N (8p~ l n N - 4p _1 - 27%) + 77^(77^ + 2p _1 - 

6p—l7rNml4 

using {7 ^ k,l,m} - {j =t k =t I = m] - {j =t k = I =t m] - 
{j ± k = m± l}-{j * k = l = m} 
where the first bracket yields 

^ sj (n) (s n - s_j (n)) (s - Sj (o)) (s o - s_j (o)) 

= ^l^nll^o I ~~ SnSoTn-o ~ SnSoTn+o-1 + P 2 l-^n I P M^ol 

_ +p- 2 ^ 2 |+p- 2 |s 2 |-p- 3 _ 

= (7jv - 2p 1 (7 JV + 3p 1 n N U N - n N U N - p 1 nj l 

= lift — 2p~ 1 U N + p~ l n N (U N - n N ) + 2p~ 1 n N U N - n N U N 

= Ujj + 2p~ 1 U N (n N - 1) + p _1 7r w (I/ w - n N ) - n N U N 



together 

U]i + 2p- 1 U N (n N + p~ l n N (U N - n N ) - n N U N 

-n 2 n N (U N - n N ) 

-p _1 (l - 2n N ){U N -n N ) 

-p _1 (l - 2n N ){U N -n N ) 

-p _1 7r w (f/ w -n N ) 
= Ufi + {-tt' 2 ti n - 2p~ 1 (l - 2tt jv ))((7 jv - n N ) + 2p- 1 U N (n N - 1) 

— n N U N 

= Vn + (U n - n N ){6p~ 1 n N -n N - 2p _1 ) + 2p~ 1 U N {Ti N - 1) 

— n N U N 

= (7jv + U N (8p~ 1 n N - 4p _1 - 27TJV) + n N (n N + 2p _1 - 6p _1 7r w ) 



Check: m k = m\ -> £"(p^) = mi(7jv 
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Orthogonal Evolution and Anticipation 



The variance is 

Var(p N ) = E(pj,) - (n N (m 2 - ml) + mlU N ) 2 

= E(pl) - n\m\ + 2n N (n N - U N )m 2 m 2 - (n N - U N ) 2 rr4 
= p~ 1 n 2 ] m 4 + 4p _1 7r JV (f/ JV - 7r w )m 1 m 3 + (n'tU^ + 2p~ 1 n N (l - n N ) - n^)m 2 
+ ((2n' 2 n N + 4p~ 1 (l - 2n N ) - 2n N )(U N - n N ))mlm 2 
+ {Ufj + U N (8p~ l n N - 4p _1 - 2n N ) + n N (n N + 2p _1 - 6p~ l n N ) 
-nl + 2n N U N - Ufj)m\ 
= p~ 1 n 2 ] m 4 + 4p _1 7r JV (f/ JV - n N )m 1 m 3 + p~ l n N {2 - 3n N )m 2 
+4p _1 (l - 3n N )(U N - n N )m\m 2 + 2p~ 1 (2U N (2n N - 1) + n N (l - 3n N ))ml 



One easily verifies 



VardpJ { 



= (m n = m\) 

6 0(p _1 ) (otherwise) 

irjv + tz n (U n — n N ) + n N (2 — 3n N ) 

+4(1 - 2n N ){U N - n N ) + 2(2U N (2n N -l) + n N (l- 3n N )) 
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